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ABSTRACT 

We describe in detail our characterization of the compact radio source population in 140 GHz 
Bolocam observations of a set of 45 massive galaxy clusters. We use a combination of 1.4 and 30 GHz 
data to select a total of 28 probable cluster member radio galaxies and also to predict their 140 GHz 
flux densities. All of these galaxies are steep-spectrum radio sources and they arc found preferentially 
in the cool-core clusters within our sample. In particular, 11 of the 12 brightest cluster member 
radio sources are associated with cool-core systems. Although none of the individual galaxies are 
robustly detected in the Bolocam data, the ensemble-average flux density at 140 GHz is consistent 
with, but slightly lower than, the extrapolation from lower frequencies assuming a constant spectral 
index. Specifically, we find a multiplicative factor of 0.85 ± 0.16 between the flux densities observed 
at 140 GHz and those predicted from a power-law extrapolation. In addition, our data indicate an 
intrinsic scatter of ~ 30% around the power-law extrapolated flux densities at 140 GHz, although our 
data do not tightly constrain this scatter. For our cluster sample, which is composed of high mass and 
moderate redshift systems, we find that the maximum fractional change in the Sunyaev-Zel'dovich 
signal integrated over any single cluster due to the presence of these radio sources is ~ 20%, and only 
~ 1/4 of the clusters show a fractional change of more than 1%. The amount of contamination is 
strongly dependent on cluster morphology, and nearly all of the clusters with > 1% contamination 
are cool-core systems. This result indicates that radio contamination is not significant compared to 
current noise levels in 140 GHz images of massive clusters and is in good agreement with the level of 
radio contamination found in previous results based on lower frequency data or simulations. 

Subject headings: cosmology: observation — galaxies: clusters: general — radio continuum: galaxies 



1. INTRODUCTION 

Recently, large Sunyaev-Zel'dovich (SZ) effect sur- 
veys from the Atacama Cosmology Telescope (ACT), 
the South Pole Telescope (SPT), and Planck have de- 
livere d catalogs with hundreds of massive galaxy clus- 
ters dSunvaev and Zerdovichl 119721: IVanderlinde et alj 
20101: IMarriage et all I2011at IRe ichardt et all I2012at 
Planck early results Villi 1201 j ). In the case of ACT 
and SPT, which operate from the ground with modest 
frequency coverage near 150 GHz and arcminute reso- 
lution, unresolved galaxies present a potential system- 
atic uncertainty in characterizing (or even detecting) 
galaxy clusters via the SZ effect flCoorav et alj H99« 
Massardi fc de Zottil [2001 : ICoble et alj 120071: iLin etaf] 
20091 : iSehgal et all I2010T ). Both instruments have ac- 
curately meas ured the blank sky population of bright 
radio sources (|Vieira et al.ll2010l IMarriage et al.ll2011b|) 
and the anisotropy power spectrum from both ra- 
dio and dusty submillimetcr sources (|Das et all I20TTI : 
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IReichardt et al.|[2012bf h However, there are few observa- 
tional constraints on the 150 GHz emission from mem- 
ber galaxies in massive clusters, even though clusters are 
known to host a significant number of radio galaxies. For 
example, within 0.5 arcmin of the cluster center there are 
~ 30 times more radio sources compared to blank sky 
and from 0.5 to ~ 2 arcmin of the cluster center there 
are ~ 3 times more radio sources compare d to blank 
sky (jCoble et all 120071 : iMuchovei et all 120101) . Unfortu- 
nately, the study of extragalactic radio sources at mm 
wavelengths is relatively immature, and it is therefore 
non-trivial to properly account for these radio sources in 
150 GHz SZ data. Even worse, the bulk of the 150 GHz 
observational data is limited to very bright flat-spectrum 
sources and is not particularly relevant to the dimmer 
steep-spectrum sources typical of central cluster galaxies 
(ICoble et al J 120071 ISadTer et alJl200^ iVieira et aijl2010t 
Marriage et all l2011bl: iGold et al.1 120111 : iSajina et aTl 
20lH\Planck early results XIHll201ll) . 

Fortunately, the low frequency properties of radio 
sources are well understood. First, radio sources are 
generally grouped into two families: flat-spectrum and 
steep-spectrum. Steep spectrum sources have a ~ 1 GHz 
spectral index of a < —0.5 and flat spectrum sources 
have a ~ 1 GHz spectral index of a > —0. 5, where a de- 
scrib es a source spectrum of the form v a (|de Zotti et al.l 
2010). Physically, a flat spectrum is typically due to 
synchrotron self-absorption (i.e., the medium is optically 
thick), which usually means the jet is oriented towards 
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the viewer. In contrast, when the jet is orthogonal to the 
viewing angle, the optically th in lobes are visible, w hich 
tend to have a steep spectrum (|de Zotti et al.ll2~010D . For 
both families of radio sources, the spectral index tends 
to decrease at higher frequencies due to electron aging 
(rclativistic energy loss) and/or the med ium becoming 
less optically thick at h igher frequencies (|de Zotti et al.l 
2010; Kcll ermannl 1 1 96(t ) . This steepening has been ob- 
served in both flat spe c trum s o urces (e.g.. iSadler et al.l 
(120081): [Marriage et al.1 (l2011bl): IMassardi et all (120101) : 
\Planck early results XIIJ (|2011[ ) HSaiina et al.l (2011 1 
and steep s pectr um sources (e.g.. iRicci et al.1 ( 2006 ): 
iTucci etahl (|2011l )). altho ugh some measur e ments do not 
find su ch steepening (e.g., iLin et all (|2009f ): iVieira et al.l 
(|2010f )). In addition, steep spectrum sources are gener- 
ally associated with powerful elliptical and SO galaxies 
(i.e., the BCG of a massive cluster), and nearly all of 
the radio sourc es associate d with clusters have a ste ep 
spectrum (e.g., ITucci et al.l (|20ll ; ICoble et al.l (|200l V 
Finally, steep spectrum sources tend to have little or 
no temporal variability (ITucci et al.l 20111: iTingav et al.l 
2003; Sa dler et al . 2006; B olton et al.ll200f 



This wealth 

of knowledge has been exploited in a range of calculations 
and simulations to estimate the amount of radio contain- 



ination in 150 GHz SZ measurements I 


Lin et al.l 120091: 


iSehgal et al.ll2010t lAndersson et al.ll2011 


). These results 



indicate that radio contamination in massive clusters at 
150 GHz should be relatively minor, but they have not 
yet been systematically verified via observational data. 

2. DATA REDUCTION 
2.1. Bolocam 

During the period from 2006 to 2012 we observed a 
set of 45 massive galaxy clusters usi ng Bolocam at the 
Caltech Submillimeter Observat o ry dGlenn et al.l 1X9981; 
Haig et aTl 120041 ; iCzakon et al.l lin prepl : iSavers et al.l 



in prepf T These clusters have a median mass of M500 — 
9 x 10 14 M and masses as low as M 500 ~ 3 x 10 14 M Q , 
which is similar to the SPT and ACT ma ss limits 
(jMarriage et alj|2011al iReichardt et~aT]|2012a|) FI All of 
the masses were computed using Ch andra X-ray data ac - 
cording to the methods described in lMantz et al.l (|2010f ). 
The median redshift of our sample is z = 0.42, and the 
median SZ S/N is 12. A complete descrip tion of these 
observations and data are g i ven elsewhere (ISavers et al.l 
120111 : ICzakon et al.l lin prepl ; ISavers et al.l lin prepf ). and 
we briefly summarize the relevant details below. For all 
of these observations, Bolocam was configured to oper- 
ate at 140 GHzE3 with a 58 arcsec full-width at half- 
maximum (FWHM) point-spread function (PSF). All of 
the cluster images are 14 x 14 arcmin squares centered 
on the cluster. The centers of the images generally have 

9 M500 is the mass enclosed within a sphere with an average 
density of 500 times the background density. 

10 The SZ emission-weighted band center of Bolocam is 
140.0 GHz. The effective band center for typical radio source spec- 
tra is slightly lower (138.7 — 139.2 GHz for the radio sources with 
— 1.5 < a < —0.5). This difference in effective band center pro- 
duces less than a 1% change in the extrapolated flux density from 
30 GHz, and is therefore not included in our analysis. Further- 
more, we note that our data are calibrated against objects with 
approximately thermal spectra, and both the source spectrum and 
the Bolo cam spectral respon se were fully accounted for in the cal- 
ibration iSavers et al,|l2012b] ). 



noise RMSs of 0.7 — 1.5 mJy/beam, (15 — 30 ^Kcmb- 
arcmin) increasing to approximately twice that value at 
the lower coverage edges of the images. T his depth is sim- 
ilar to the SPT and ACT surv ey depths (jMarriage et al.l 
l2011atlReichardt et al.|[2Tjl2al) . 

In order to subtract atmospheric fluctuations from 
these data, we remove the time-instantaneous average 
signal over the field of view and then high pass filter the 
datastreams at a characteristic frequency of 250 mHz. 
This process also removes some astronomical signal from 
the data, and we quantify this via simulation to obtain 
a signal transfer function. For the results presented in 
this manuscript, the cluster signal transfer funct i on wa s 
computed as described in detail in ISavers et al.l (|2011l) . 
We did not explicitly compute transfer functions for the 
point sources. Instead, we processed point source mod- 
els centered on the galaxy positions determined from low 
frequency data through our reduction pipeline and then 
directly fit the resulting filtered image to our data. In all 
cases the point source models were normalized to an am- 
plitude of 1 mJy. We note that the effective filtering of 
the point source model is independent of the amplitude 
of the model for sources < 100 Jy, which is at least 3 or- 
ders of magnitude brighter than any of the radio sources 
in our sample. 

2.2. OVRO/BIMA and SZA 

Most of our radio source flux densities at ~ 30 GHz 
were obtained from previously published results from 
the Owens Valley Radio Observatory (OVRO), the 
Bcrkclcy-Illinois-Maryland Association (BI MA) array, 
and the Sunyaev-Zerdovich Array (SZA) (jCoble et al.l 
l2007tlBonamente et al.ll2012l ). However, we reduced data 
collected by these facilities for an additional 13 radio 
sources (3 from OVRO/BIMA and 10 from SZA). Our 
redu ction of these data follo wed the methods describ ed in 
e.g.. ICoble et al.l (|2007t) and lMuchovei et al.l (|2007t) . and 
we refer the reader to those manuscripts for additional 
details. To obtain point source flux densities we simulta- 
neously fit for the cluster SZ signal and all of the point 
sources. We modeled the clu ster SZ emissi o n usin g the 
best-fit gNFW profiles from lArnaud et all (I2010D . and 
corrections for the SZA primary beam were accounted 
for in both the cluster and point source fits. 

3. CLUSTER-MEMBER RADIO GALAXIES 

3.1. Central Radio Galaxies 

There are a significant number of bright radio galax- 
ies nea r the centers of massive clusters (e.g.. ICoble et al.1 
(|2007f) find that ~ 97% of radio sources within 0.5 ar- 
cmin of the cluster center are cluster members). Un- 
fortunately, these sources are difficult to detect in our 
arcmin-resolution SZ data due to the large degeneracy 
between the inner slope of the SZ profile and the radio 
source flux density. Consequently, we have searched the 
1.4 GH z NVSS catalog wit hin 1 arcmin of each cluster 
center (jCondon et al.l fl998) in order to select probable 
cluster member galaxies E3 We chose the NVSS survey 

11 Throughout this paper we refer to the NVSS selected sources 
within 1 arc min as cluster mem ber galaxies. However, based on 
the results of Coble et al. (2007) and the blank-sky source density 
of the NVSS catalog, we expect > 5% of the sources to be field 
galaxies unassocated with the clusters. 
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Fig. 1. — Left: the number of NVSS sources as a function of radius from the cluster centers. We show the number of NVSS sources for 
the full sample of 45 clusters, the 16 disturbed clusters, and the 17 cool-core clusters. Note the significant number of NVSS sources at 
the very center of cool-core clusters. For reference, we also show the blank-sky average number of sources for the NVSS survey given the 
angular size of each bin. Right: cluster masses for both the full sample of 45 objects and for the su bset of 25 ob j ects w ith central radio 
sources. The masses were computed from Chandra X-ray data according to the methods described in [Mantz ct al. (201JJ). 

best-fit SZ and point-source templates into 100 noise re- 
alizations and repeated the fit in order to characterize our 
uncertainty on the point-source flux density. For most of 
our data the raw point-source sensitivity is ~ 1 mJy, 
but our uncertainties on the flux densities of these point 
sources are ~ 1 — 3 mJy, with the degradation due to de- 
generacies between the SZ and point source templates. 
The sensitivity degradation depends strongly on the sep- 
aration of the SZ and point-source centroids, along with 
the angular size of the SZ template. For clusters with 
multiple radio sources, we fit for each source separately. 
The best-fit flux densities are given in Table [2] 

We then derived an independent estimate for the 
140 GHz flux densities of the sources by extrapolating 
the low frequency data. We have primarily combined 
1.4 GHz NVSS data with ~ 30 GHz measurements from 
SZA and/or OVRO/BIMA, which exist for 24/28 sources 
(IBonamente et all |20T1 fCoble et alJl2007h FI Three of 
the sources without 30 GHz data have measurements 
at either 4.85 GHz and/or frequencies below 1.4 GHz, 
which have been used to constrain t he spectral energy 
distributions (SEDs) of those sources (iCohen et al.ll2007t 
ILaree et al.l[l98l IGriffith et all 11991 . We fit an SED 
of the form v ai - 4 /3o, and we find all of the galaxies have 
a i.4/30 < —0.50 (see Figurel2"j) F^l Although our source se- 
lection using 1.4 GHz data favors steep spectrum sources, 
we note that the ~ 30 GHz data do not contain any 
point source detections without an NVSS counterpart. 
The ~ 30 GHz data are generally sensitive to sources 
with flux densities larger than 0.5 — 1.0 mJy, and would 
therefore have detected bright flat spectrum sources if 



because it has approximately uniform depth over the full 
sky above dec = —40, and is complete abov e a flux den- 
sity of 2.5 mJy. Motivated by the results of ICoble et al.l 
(|2007| ) , the 1 arcmin search radius was chosen to be large 
enough to contain the bulk of the bright cluster member 
radio galaxies without containing a significant number of 
foreground and background galaxies that are not associ- 
ated with the cluster (see Figure [T]) . Over our sample of 
45 clusters, we find 31 NVSS sources spread between 26 
clusters (three clusters contain two sources and one clus- 
ter contains three sources, see Table [T j) . In agreement 
with previous results (e.g.. iMittal et a l. (2009)), most of 
our cool-core clusters have at least one compact NVSS 
source (15/17), while less than half of the disturbed clus- 
ters in our sample have a compact NVSS source (6/16) 03 
We note that two systems are both cool-core and dis- 
turbed, so the number of non-cool-core disturbed systems 
with compact NVSS sources is 4/14. 

For each of these NVSS selected sources we have mea- 
sured the 140 GHz flux density as follows. First, we 
jointly fit a point source template centered on the NVSS 
coordinates, along with a template for the SZ emission, 
to our data. Our fitting procedure is described in de- 
tail in ISavers et al.l (| 2 1 If) , and we refer the reader to 
that manuscript for additional details. The normaliza- 
tion of the point-source template gives the flux density 
of the source (in mJy). For the SZ template, we used 
the morphology-dep endent best-fit pressure profiles from 
lArnaud et al.l (|20Tol ) (e.g., we used the best-fit cool-core 
template for the cool-core clusters in our sample), allow- 
ing the normalization of the template (i.e., po) to be a 
free parameter. These profiles were used because all of 
the point sources we fit are in the central regions of the 
cluster, where the shape of the pressure profiles varies sig- 
nificantly with morphological type. We then inserted the 

12 Our classification of clusters as cool-core or distu rbed is based 
on X-ray measurements and is described in detail in Saycrs ct al. 
I|in predr i . Briefly, the projected luminosity ratio is used to define 
cool-corc systems (Mantz 2009; Mantz ct al. 2010; Bohringcr ct al. 
120101) an d the centroid shift parameter is used to define disturbed 
syste ms UMaughan et a l. 2008; Pratt et all 120091 ; IMaughan et all 
120121) . 



13 3/31 of the NVSS sources are extended and appear to have 
very steep spectral indexes. Extrapolating to 140 GHz, we expect 
negligible flux density from these three sources and therefore have 
not attempted to fit to them in our data. 

14 For our SED fits we assume a single observing frequency 
for each dataset independent of source spectrum. This results in a 
negligible bias in our results due to the small fractional bandwidths 
and large separations in band centers among our datasets (1.4, 30, 
and 140 GHz). Furthermore, we note that there are some system- 
atic differences in calibration between the datasets (e.g., using prc- 
or post-WMAP planetary models), but these differences are below 
the flux calibration uncertainties quoted for each dataset. 
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TABLE 1 

NVSS SOURCES WITHIN 1 ARCMIN OF THE CLUSTER CENTER 



cluster 


rcdshift 


morphology 
cool-core disturbed 


RA (J2000) 


dec (J2000) 


NVSS sources 

distance (arcsec) 


flux density (mjy) 


MAOo J0018.5 


0.54 








no sources 




MACb J0025.4 


0.58 




t \!\ of oo f\ 

00:25:32.0 


TO OO OO "1 

— 12:23:09.1 


39 


O O T 1 "1 O 

28.7 ± 1.3 


ZWtL UU24 


0.39 


/ 


00:26:35.5 


+ 17:09:31.7 


10 


2.6 ± 0.4 


Abcll 209 


0.21 




01:31:52.7 


10 1 ) / ' r" n 

— 13:36:59.6 


13 


18.4 ± 1.0 


CL J0152.7 


0.83 


/ 






no sources 




A 1 11 OC7 

Abell 2o7 


0.23 


/ 






no sources 




A i n oTn 

Abcll 3 i 


0.38 


/ 


no or* c rr o 

02:39:55.3 


m 1 a 1 1 

— Ul:o4:lo.o 


37 


10.5 ± 1.0 


Abcll ooo 


0.19 




no. a . no a 
02:4o:0o.4 


n q . q 1 . a q n 
— U3:31:43.9 


3 


.inn _l_ 1 1 
4U.y ± 1.3 


IVlAL/b JU2o ( .1 


0.50 








no sources 




i\ .f A no irnon c; 

JVlAGb JUozy.u 


0.45 


/ / 


no on a 1 t 

03:29:41.7 


no 11 coo 

—02:11:52.2 


7 


6.9 ± 0.6 


MAGb JU410.1 


0.42 


/ 






no sources 




JVlAGb JU417.5 


0.44 


/ / 


f\ A 1 T O A n 

04:17:34.9 


— 11:54:34.2 


12 


on 1 1 t 

29.8 ± 1.7 


i\ .f A /^i o in ,1 on c; 


0.40 


/ 


c\a . on . n 
04:29:30.0 


— U2:o3:U0.4 





1 QO O _l_ A O 

loo.o ± 4.2 


ii j Ann Tn/ici n 

MAOo JU451.9 


0.43 


/ 






no sources 




|\» A rin \~f~\AVi A 1 

MAub JU454.1 


0.55 


/ 






no sources 




i\ if A no Tr\r" a t 

MAOa J0647.7 


0.59 








no sources 




MACS JU717.5 


0.55 


/ 


n^7. 1 *7 . Q A 1 

U7:l ( :34.1 


+0 f : 45:01. 3 


31 


nn n _i_ q *7 

yu.y ± 0. ( 








n'y 1 7 op •? 

07:17:35.7 


i >-7 ^ c on 

+37:45:39.2 


47 


1 no rr I o *7 

1U2.5 ± 3.7 


MACb JU744.5 


0.69 


/ 






no sources 




Abcll oil 


0.29 








no sources 




A i ii rr\ r 7 

Abcll 697 


0.28 








no sources 




JVlAGb JU911.2 


0.50 








no sources 




A i ii nc o 

Abcll 9o3 


0.21 








no sources 




Mb !Uo4 


0.83 


•/ 


in rr rr n 

10:50:59.0 


no "7 

—03:37:26.8 


18 


14.1 ± 0.9 


MAOb Jlllo.o 


0.36 


/ 


11:15:51.8 


+01:29:55.5 


2 


16.2 ± 1.0 


MAub J1149.5 


0.54 








no sources 




Abcll 14zo 


0.21 




11:57:16.8 


1 OO Off A A C 

+33:36:44.6 


9 


00 1 n n 

33.3 ± o.y 


MAub J lzUb.z 


0.44 




1 n/? 101 

12:00:12.1 


no a no rr. 

— U8:4o:(J2.o 
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t cr\ n 1 c 


pT T i one 


0.89 




1 0. o^; . c q 
lz:zO:oo.o 


1 Q Q . Q O . A A Q 

+oo:o2:44.o 
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4.3 dz 0.5 


MAub J1311.U 


0.49 


/ 






no sources 




lv t A no n r 

MACb J 1347.5 


0.45 


/ 


I /i ^ 00 I 
13:47:30.1 


1 I (i r n 

— 11:45:30.2 


23 


17-,^- | OO 

17.7 ± 3.2 








10 A *7 on "7 

13:47:30.7 


1 1 /i rr no c 

— 11:45:08.6 


2 


45.9 ± 1.5 


Al 11 "lOOP" 

Abcll 1835 


0.25 


/ 


"1 A OI OO "1 

14:01:02.1 


1 no pf /i "i n 

+02:52:41.0 


2 


on . ) 1 / • 

39.3 ± 1.0 


ATA t~AO T 1 A OO O 

MACb J 1423. 


0.55 


/ 


1 a 00 a <~7 n 

14:23:47.9 


1 a r\ a onn 

+24:(J4:v>y.y 


3 


8.0 ± 1.1 


TNyTAr^Q 11 toil n 


0.36 




1K.QO.KQ O 

lo:32:53.o 


1 Qn.on.Kn 
+oU:2U:oy.o 


1 


OO 1 n 
22.5 ± U.o 


Abcll 2204 


0.15 


/ 


16:32:46.9 


1 n cr /i ,1 n 

+05:34:34.9 


5 


f n 1 k 

69.3 ± 2.5 


Abcll 22ly 


0.23 




1 c /in 1 er n 

16:40:15.0 


1 /If -1 O OO >"7 

+46:42:28.7 


55 


6.1 ± 0.5 








1 c A n 1 

16:40:21.8 


1 AC -IO A *7 O 

+46:42:47.8 


24 


oon 1 I 

239.1 ± 8.3 








16:40:23.8 


+46:41:47.3 


56 


7.9 + 1.0 


MACS J1720.3 


0.39 


/ 


17:20:15.3 


+35:35:26.3 


59 


9.6 + 0.5 








17:20:16.8 


+35:36:28.4 


6 


18.0 + 1.0 


Abell 2261 


0.22 


/ 


17:22:27.7 


+32:07:57.8 


8 


5.3 + 0.5 


MACS J1931.8 


0.35 


/ 


19:31:49.9 


-26:35:13.4 


40 


216.5 ± 6.9 


MS 2053 


0.58 


/ 






no sources 




MACS J2129.4 


0.59 


/ 






no sources 




RX J2129.6 


0.24 


/ 


21:29:40.0 


+00:05:23.0 


7 


25.4 + 1.2 


MS 2137 


0.31 


/ 


21:40:15.0 


-23:39:39.5 


2 


3.8 + 0.5 


MACS J2211.7 


0.40 


/ 






no sources 




MACS J2214.9 


0.50 


/ 


22:14:56.5 


-14:00:55.7 


46 


5.6 + 0.6 


Abell S1063" 


0.35 








no sources 





Note. — NVSS sources within 1 arcmin of the cluster center for the 45 clusters in our sample. The columns give the cluster name, the 
cluster rcdshift, the X-ray morphology, the RA/dec of the NVSS source, the distance from the cluster center to the NVSS source, and the 
NVSS flux density at 1.4 GHz. The superscript a denot es that Abell S1063 is below the decl i nation range of the NVSS survey (—40 deg), 
and we therefore searched the SUMSS catalog instead IIBock et al.| [T999: Mauch ct al. 2003: IMurphv et al .1120071 '). The SUMSS survey is 
similar in depth to the NVSS survey, although at lower frequency. We found no SUMSS sources within 1 arcmin of the center of Abcll 
S1063. 



they exist in the central regions of these clusters. Since 
all of the bright central radio sources are steep spectrum 
sour ces, which tend to have little or no temporal variabil- 
ity dTingav et al.ll2003t iSadler et all 120061: iBolton et al.l 
|2006J), the asynchronous NVSS/SZA/OVRO/BIMA ob- 
servations should provide an accurate characterization of 
the source SED. We used the values of ai.4/30 to predict 
the flux density at 140 GHz, and these extrapolated val- 
ues are given in Table [TJ 

None of the NVSS sources near cluster centers are de- 
tected at a significant level in our 140 GHz images. Con- 
sequently, we have characterized the average point-source 



signal measured in our data by comparing our measured 
flux density to the predicted flux density based on the 
extrapolation from lower frequency measurements. To 
avoid correlations and biases in this average measure- 
ment, we have discarded the data from clusters with mul- 
tiple NVSS sources, those with no 30 GHz data, and also 
Abell 2261 due to the fact that the SZ template does 
not provide a good description of those data. For the 17 
remaining clusters we formed the quantity 



SboIo — S, 



cxtrap 



S/Bolo 



S I cxtrap 
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TABLE 2 
Central Radio Sources 



cluster 


NVSS 


OVRO/BIMA 


SZA 


a 1.4/30 


extrap. 


Bolocam 




1.4 GHz 


28.5 GHz 


30.9 GHz 


140 GHz 


140 GHz 




mJy 


mJy 


mJy 




mJy 


mJy 


MACo JUU25.4 


28.7 ± 1.3 


n o o I n 1 tc 

0.33 ± 0.17 




— 1.38 


U.Uo rb U.Uo 


—2.3 ± 2.0 


ZiWL'ij UU24 


2.6 ± 0.4 




— 0.2 dr. 0.2 


— 1.09 


n m _L n no 
U.U1 ± U.U2 


0.6 ± 5.3 


Abcii zuy 


18.4 ± 1.0 




1.1 rb 0.2 


—0.92 


n _i_ n 117 
U.zo ± U.U-r 


— 1.2 ± 2.1 


AT n Q r7(-i 

Abcii o7U 


10.5 ± 1.0 


U.72 ± U.U8 


0.8 dz 0.2 


—0.89 


n 1 1 n r\i 

U.lo ± U.Uo 


—3.4 ±1.8 


A QOQ 

Abcii ooo 


40.9 ± 1.3 


4.4U rb U.4 ( 


A Q _L n o& 
4.o ± U.2 


—0.73 


i /in _i_ n 1 n 
1.4U ± U.1U 


5.4 ± 2.1 




6.9 ± 0.6 




0.3 dz 0.4 


— 1.04 d 


n 1 n I n 10 

U.1U dr. U.lz 


3.9 ±3.1 


TV/TAf^Q 70/11*7 K 


on o _l_ 1 *7 
29.8 ± 1. / 






— 1.11 


n 1 _L n no 
U.lo ± U.Uo 


2.3 ± 2.1 


TV/TA/^Q Tri/ion £; 
MAU& JU4zy.u 


IOC O _1_ A O 

138.8 ± 4.2 




i o o _L n oi> 
lo.z ± U.Z 


—0.66 


u.bo ± U.lo 


9.8 ± 3.5 


TV JQ 1 n C /l 

Mb !Uo4 


14.1 ± 0.9 


n c\a _i_ n ncQ 
u.y4 ± U.Ub 




—0.90 


n 00 _L n no 
U.zz rb U.Uz 


1.6 ± 1.7 


1\/T A /^G 7111c; Q 
ivlAL/O Jlllo.o 


1 a o _l_ 1 n 




1 'i i n 4b 
1.4 rb U.4 


— U.ol 


n /io 1 n 10 
U.4Z rt U.lo 


7 _i_ 1 
— z. ( ± i.y 


Abcii 14zo 


ooo I n n 

oo.o ± u.y 




0.6 zb 0.2 


— 1.31 


U.U9 ± U.U4 


—3.5 ± 2.8 


Tv/TAr^Q 71 one o 
MAU& JlzUb.Z 


ibu.y ± O.O 






— 1.38 d 


n 0*7 _L n no 
U.2 / ± U.Uo 


—5.4 ±1.9 


L/L/ Jlzzb 


4.3 rb 0.5 




0.3 dr. 0.2 


—0.92 


n nn _i_ n no 

u.uy ± u.uo 


1.3 ± 3.3 


MAtb Jlo4/.D 


45. y ± 1.0 


1 n q q _L n 1 t: a 
1U.OO ± U.15 


t j_ n oi> 
0. / ± U.2 


—0.51 


/i /i _l_ n 1 1 e 
4.4Z rb U.ll 


— 4.4 ±1.8 


Abcii loob 


on Q _L 1 d 

oy.o ± l.o 


z.y / ± U.Ub 


z.y ± U.o 


—0.86 


n *7*x _L n no 
U. 10 ± U.Uo 


— 1.1 ± 1.4 


ATA /~1 U irin o 

MAub J14zo.o 


8.0 rb 1.1 


1.49 rb U.lz 


n 1 n 06 

2.U ± U.2 


—0.51 


rv r 7fz 1 rv 1 rv 

U.75 ± U.1U 


2.1 ± 2.8 


T\I A PC 11 coo n 


OO Q I n Q 

zz.o rt U.o 


q oc; _l n 1 qq 
o.zo rt U. lo 


Q O J_ A 0& 

o.z rb U.Z 


— U.O 1 


1 17 x n n7 
1.1 ( rt U.U / 


r 1 O r? 

z.o ± z. ( 


Al TT oon (I 

Abcll 2204 


69.3 ± 2.5 


8.79 ± U.lo 


7.0 ± 0.2 


—0.70 


2.75 ± 0.07 


1.1 ± 1.4 


Abell 2219 


6.1 ±0.5 




0.6 ±0.2 


-0.74 


0.21 ±0.09 


4.9 ±2.6 




239.1 ±8.3 


14.87 ±0.17 a 


14.2 ±0.2 


-0.92 


3.42 ± 0.07 


-4.9 ±3.3 




7.9 ± 1.0 




0.1 ±0.1 


-1.40 


0.02 ±0.02 


-7.3 ±2.4 


MACS J1720.3 


9.6 ±0.5 




0.6 ±0.1 


-0.90 


0.15 ±0.04 


5.3 ±1.6 




18.0 ± 1.0 




1.8 ± 0.4 6 


-0.76 


0.58 ±0.21 


-5.1 ±2.4 


Abell 2261 


5.3 ±0.5 


0.20 ±0.30 




-1.10 


0.05 ±0.07 


-4.8 ±1.3 


MACS J1931.8 


216.5 ±6.9 






-0.72 d 


8.81 ±4.56 


-3.0 ±2.2 


RX J2129.6 


25.4 ± 1.2 


2.33±0.10 a 


2.6 ± 0.2 6 


-0.77 


0.70 ± 0.04 


0.1 ±2.9 


MS 2137 


3.8 ±0.5 






N/A 


0.06-f 


-7.5 ±5.1 


MACS J2214.9 


5.6 ±0.6 




0.4 ±0.3 


-0.89 


0.14 ±0.13 


-2.8 ±2.3 



Note. — Compact central radio sources in our sample of 45 clusters. The columns give the cluster name, the NVSS flux density at 
1.4 GHz, the OVRO/BIMA flux density at 28.5 GHz, the SZA flux density at 30.9 GHz, the spectral index determined from the 1.4 
and ~ 30 GHz data, the predicted flux density at 140 GHz based on that spectral index, and the measured flux density in our 140 GHz 
image at the location of the source, which is in some cases negative due to noise fluctuations. All of the uncertainties are lcr values; the 
Bolocam uncertainties are approximately Gaussian and in clude degeneracies between the point source and the cluster SZ signa l (see text). 
Superscripts denote the followi ng: a) measurements from [Coble et aLl 1120071 ). b) measurements from Bonamcntc et al. (2012), c) sources 
published in Coble et al. (2007), but refit in our analysis (the source associated with MACS J0025. 4 was published in B1950 coordinates 
in Coble ct al. (2007)), d) spectral index based on fits with 4.85 GHz and/or lower frequency data (Cohen et al. 2007; Large et al. 1981; 
IGriffith ct al .119941) . e) source detected by two groups at mm/submm frequencies , and the best-fit 140 GHz flux density from the combination 
of those data and the radio results is 4.4 ± 0.3 mJy ( Pointccoutcau et al. 2001) and 5.5 ± 0.6 mjy (Komatsu ct al. 1999), consistent with 
our extrapolated estimate, /) only 1.4 GHz data was available for this source, so it was extrapolated based on the median value of O1.4/30 
for our sample (-0.89). We note that some of the extrapolated flux densities arc consistent with 0. We quote symmetric uncertainties about 
these best fit values, even though the true probability density functions are asymmetric and exclude unphysical negative values. 



where SboIo is the best-fit flux density in our Bolocam 
image, S^xtrap is the extrapolated flux density based on 
the lower frequency measurements assuming a constant 
spectral index, CT5/B0I0 is the uncertainty on the Bolo- 
cam measurement, and Ps/cxti-np is the uncertainty on 
the extrapolated flux density. 

If our extrapolation of the spectral index found at lower 
frequencies perfectly describes our data then the values of 
8 will have a mean of and standard deviation of 1. We 
find a mean S = —0.26 ± 0.28, indicating that our data 
show no statistically significant evidence for a change 
in the spectral index above 30 GHz, although our data 
are consistent with a slight steepening of the spectral in- 
dex (a steepening is consi stent with previou s results and 
theoretical expectations (|Tucci et al.l 120111 : iRicci et aLl 
I2006L |Kcllermaiin 1966])). Furthermore, we note that 
a Kolmogorov-Smirnov test shows that our values of S 
are consistent with a Gaussian distribution with a prob- 
ability to exceed of 0.60. Our constraint on 8 provides 
an estimate for the average difference between the mea- 



sured and extrapolated 140 GHz flux densities relative 
to the uncertainty on the average flux density. We there- 
fore convert these constraints on S to constraints on the 
extrapolated flux densities via the noise-weighted mean 
140 GHz flux density for these sources (0.93 ±0.52 mJy). 
Based on this mean flux density and our mean measured 
8, we find an ensemble- average multiplicative factor of 
0.85 ±0.16 between our 140 GHz measured flux densities 
and the extrapolated flux densities based on a constant 
spectral index. Furthermore, we find that the standard 
deviation of the 17 values of 8 is equal to 1.14 ± 0.18. If 
there were no intrinsic scatter in the extrapolation from 
lower frequencies then we would expect a standard devi- 
ation of 1. Our data are consistent with this value, al- 
though they indicate a best-fit intrinsic scatter of ~ 30% 
in the extrapolated 140 GHz flux density. 

Since the extrapolation of the radio galaxy flux densi- 
ties from lower frequencies provides a good description 
of our 140 GHz data, we subtracted each of the radio 
galaxies using these extrapolated flux densities multi- 
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TABLE 3 
Compact Central Radio Sources 



cluster 


subtracted PS 


SZ S/N (raw) 


8 C (raw) 


SZ S/N (PS-sub) 


0c (PS-sub) 


A S/N 




mjy 




arcmin 




arcmin 


percent 


TV T A o Tnrvo rr a 

MALb JUUz5.4 


0.05 


12.31 


0.50 


12.33 


0.50 


0.2 


ZWCL 0024 


0.01 


3.26 


3.00 


3.26 


3.00 


0.0 


Abcll 21)9 


0.24 


13.93 


1.75 


13.99 


1.50 


0.4 


a i n <~)Tr\ 

Abcll o7U 


0.15 


12.80 


0.50 


12.88 


0.50 


0.6 


AT 1 1 onoo, 

Abcll ooo 


1.19 


9.33 


2.00 


9.56 


1.75 


2.5 


tv /t a *~*c Trvoon c 


0.09 


12.09 


0.25 


12.13 


0.25 


0.3 


TV /T A /~>0 TrV/ll 1 ? c 

MAUb JU417.5 


0.15 


22.67 


1.00 


22.71 


1.00 


0.2 


MACS J0429.6 


5.66 


7.34 


1.50 


8.78 


1.00 


19.6 


MS 1054 


0.19 


17.38 


0.50 


17.53 


0.50 


0.8 


tv iT A nc t 1 i i r o 

MACS J1115.8 


0.36 


10.92 


1.25 


11.03 


1.25 


1.0 


Abell 1423 


0.08 


5.80 


0.50 


5.82 


0.50 


0.3 


TV X A /""lO T "1 nA/ 1 o 

MACS J 1206.2 


0.24 


21.69 


1.00 


21.79 


1.00 


0.5 


/""it T 1 one 

J 122o 


0.08 


13.01 


0.25 


13.05 


0.25 


0.3 


MALb Jlv>4/.0 


3.76 


34.04 


0.50 


36.50 


0.25 


i .2 


A 1 11 lOOf 

Abell looo 


0.64 


15.34 


0.75 


15.70 


0.50 


2.4 


TV iT A nc T 1 A no o 

MACS J 1423.8 


0.64 


8.98 


0.50 


9.34 


0.50 


1.0 


MACS J1532.9 


0.99 


7.55 


0.75 


7.98 


0.50 


5.7 


Abell 2204 


2.34 


20.94 


1.00 


22.24 


0.75 


6.2 


Abell 2219 


2.91 
0.18 
0.02 


10.53 


1.00 


11.04 


1.00 


1.8 


MACS J1720.3 


0.49 
0.13 


10.32 


0.25 


10.61 


0.25 


2.8 


Abell 2261° 


0.04 


10.76 


0.25 


10.79 


0.25 


0.3 


MACS J1931.8 


7.49 


9.21 


1.25 


10.04 


1.00 


9.1 


RX J2129.6 


0.60 


7.76 


0.75 


7.96 


0.75 


2.6 


MS 2137 


0.05 


6.51 


0.25 


6.53 


0.25 


0.3 


MACS J2214.9" 


0.12 


12.85 


1.00 


12.86 


1.00 


0.1 



Note. — Compact central radio sources subtracted from Bolocam cluster data, assuming that the flux density extrapolated from lower 
frequencies should be multiplied by 0.85 to account for our best-fit to the spectral steepening at high frequencies. From left to right the 
columns give: the flux density of the subtracted radio source, the SZ S/N in our unsubtracted image, the filter scale corresponding to this 
peak SZ S/N in our unsubtracted image, the SZ S/N in our point-source subtracted image, the filter scale corresponding to the peak SZ S/N 
in our point-source subtracted image, and the fractional difference in SZ S/N between the unsubtracted and point-source subtracted images. 
In general, the extrapolated flux densities of the radio sources are accurate to ~ 30%, limited by the intrinsic scatter in the extrapolation. 
Consequently, the fractional change in the SZ S/N is also accurate to ~ 30%. The superscript a denotes clusters that contain point sources 
directly detected in our Bolocam images away from the cluster centers (see Section . These additional point sources were not subtracted 
for this analysis. 

radio sources is given in Table |3] In most cases the mea- 
surement uncertainty on the extrapolated flux density 
is quite small, and our overall uncertainty is therefore 
dominated by the 30% intrinsic scatter in the extrapo- 
lation. We then computed the peak SZ S/N for each 
cluster both with and without the radio source(s) sub- 
tr acted using the optim al filtering formalism described 
in I Savers et all (|2012aT i. Briefly, we construct an SZ 
template of the form S{8) oc (1 + 2 /0j?) -1 , weight this 
template by the signal transfer function of our map and 
the inverse of the map variance, convolve our map with 
the resulting filter, and search for the peak S/N. This is 
repeated for 9 C = 0.25, 0.50, .., 3.00, and the maximum 
peak S/N over all filter scales is used to describe the 
SZ S/N. This is essentially the same algorithm used by 
both AC T and SPT to search for clusters in their sur- 
vey data (IVanderlinde et al.l[2010t iMarriage et alJl20lTat 
iReichardt et al.ll2012a|) . We find the maximum fractional 
change in SZ S/N for any of the clusters is ~ 20%, and 
only 12/45 of the clusters show a fractional change in 
SZ S/N larger than 1% (see Figure [3}. Consequently, 
for the high-mass and moderate-redshift clusters in our 
sample contamination from radio emission is in general 
small compared to the noise in our data. Furthermore, 
for the brightest radio galaxies we generally find that the 
value of 6 C with the maximum SZ S/N is smaller when 
the radio source is subtracted from the map, which is 
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all clusters, <z> = 0.42 
disturbed, <z> = 0.52 
cool-core, <z> = 0.36 
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-1.2 



■1.0 



-0.8 -0.6 



a 



1 .4/30 



Fig. 2. — The value of the spectral index O1.4/30 found for the 
cluster-member galaxies in our sample using measurements at u < 
30 GHz. All of these galaxies are steep spectrum sources, and we 
find a median spectral index of —0.89. The spectral indices for the 
disturbed and cool-core sub-samples are also shown, and have been 
offset by 0.01 for ease of viewing. Compared to the full sample, the 
galaxies associated with cool-core clusters have shallower spectral 
indices on average. 

plied by 0.85 (to account for our best-fit steepening of 
the spectral index). A complete list of the subtracted 
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Fig. 3. — Top left: The fractional change in the SZ S/N as a function of the extrapolated 140 GHz flux density of the radio source. 
The maximum change is ~ 20% and a flux density > 0.5 mJy is required to produce more than a 1% change in SZ S/N. The data are 
color-coded by X-ray derived morphology, and nearly all of the brightest radio sources are associated with cool-core systems. Over half of 
the radio sources arc within 10 arcsec of the cluster center (25 — 75 kpc for the redshift range of our sample), and these sources are marked 
as filled symbols. Note that roughly half the clusters in our sample contain no radio sources given our selection criteria, and are therefore 
not included in this plot. Top right: The cumulative fraction of clusters with radio contamination below a given threshold (as quantified 
by fractional change in SZ S/N). The black, red, and blue curves represent the full sample of 45 clusters, the disturbed sub-sample, and the 
cool-core subsample. The radio contamination is severe enough to cause a > 1% change in the SZ signal for ~ 25% of the clusters in our 
sample, although the contamination is significantly worse in cool-core systems. Note that the cumulative fraction ignores sources undetected 
by NVSS, and is therefore likely to be an overestimate at A S/N < 1%. Bottom left: The fractional change in SZ S/N versus cluster mass, 
showing no evidence for a correlation and indicating that the fractional amount of radio contamination is approximately independent of 
mass for these high mass clusters. Bottom right: The X-ray projected luminosity ratio L ra t = L(R < 0.05Rsoo)/L(R < R500) versus the 



in the clusters with the largest projected luminosity ratio. We define cool-core systems as clusters with L rat > 0.17 



re preferential! 
HMantzll2009ft . 



due to the resulting cluster profile being more sharply 
peaked (recall that these radio sources partially fill in 
the SZ decrement at 140 GHz). 

Although the overall level of radio contamination in 
our sample is small, we note that 11/12 of the bright- 
est radio sources are associated with cool-core systems, 
and therefore 11/17 cool-core systems show a fractional 
change in SZ S/N larger than 1%P^I Previous studies 
have found a similar relationship be tween stron g radio 
emission and cool-core systems (e.g., iSunl (|2009f )). Fur- 
ther supporting this trend of cool-core systems hosting 

15 The median redshift of the cool-core systems in our sample 
(z = 0.36) is slightly lower than the median redshift of the full sam- 
ple (z = 0.42) and the median redshift of the disturbed systems 
(z = 0.52). Consequently, some of the difference in radio contami- 
nation for the different morphological classifications may be due to 
these differences in redshift rather than the cluster environment. 



the brightest radio galaxies, we also find a correlation 
between the strength of the radio emission and the pro- 
jected X-ray luminosity ratio, which we have used to sep- 
arate cool-core and non-cool-core systems (see Figure [3]). 
This correlation is quantified by a Spearman rank co- 
efficient of 0.45 and null hypothesis probability of 0.03. 
In addition, we find that there is no significant correla- 
tion between the fractional change in SZ S/N and cluster 
mass for our sample, quantified by a Spearman rank co- 
efficient of -0.11 and a null hypothesis probability of 0.62 
(see Figure |3|). Consequently, the fractional change in SZ 
S/N due to radio contamination is approximately inde- 
pendent of mass for these high mass systems. Since there 
is no evidence for a correlation between cluster mass and 
fractional change in SZ S/N, and since our sample ap- 
proximately spans the mass rang e of objects discovered 
by the ACT and SPT surveys ([Marriage et all 1201 lat 
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iReichardt et al.l l2012af ) , the amount of radio contami- 
nation in those surveys should be similar to what we 
have found for our sample. However, we caution that 
those surveys contain a significantly larger fraction of 
low-mass clusters, and the negative Spearman rank coef- 
ficient indicates that our data favor (at low significance) 
an on-average increase in fractional change in SZ S/N 
with decreasing cluster mass. 

Due to the way we have selected the cluster-member 
radio sources, anything below the NVSS detection limit 
(2.5 mJy at 1.4 GHz) will not be included in our anal- 
ysis. However, even if we assume the smallest magni- 
tude spectral index found for any of the radio galaxies 



in our sample (0^ 



4/30 



-0.51), a source undetected 



by NVSS would have a maximum 140 GHz flux den- 
sity of 0.25 mJy. Therefore, it is unlikely that a source 
undetected by NVSS would have more than a 1% ef- 
fect on the optimally filtered peak SZ signal we measure 
for any of the clusters in the Bolocam sample. We do 
note that, of the 21 central radio sources detected by 
ICoble et al.l (|2007). there were two sources with spec- 
tral indices smaller in magnitude than ai.4/30 = —0.5, 
and one source with a spectral index of ai.4/30 = —0.14. 
Since our ~ 30 GHz data are most sensitive to these flat 
spectrum sources, we consider the the extreme case of a 
source with a 1 A / 30 = —0.14 and a ~ 30 GHz flux density 
just below our detection limit of < 1 mJy. Such a source 
could go undetected in our ~ 30 GHz data and still have 
a 140 GHz flux density of ~ 0.5 mJy. Such a source would 
likely produce a ~ 1% change in the SZ S/N measured 
by Bolocam (see Figure [3]). However, given the rarity of 
such flat spectrum sources in the centers of clusters, cou- 
pled with the fact that the source would also need to be 
close to our detection threshold, it is unlikely that there 
are a significant number of such sources in the Bolocam 
cluster sample. 

3.2. Non-Central Radio Galaxies 

Although the number of cluster radio galaxies drops 
quickly with radius, there are likely to be some cluster 
member radio galaxies outside of our 1 arcmin diameter 
search radius. For example, the total number of radio 
sources in the central ~ 2 arcmin radius of our images is 
enhance d by a factor of ~ 3 compared to the bl ank-sky 
average ([Coble et al.ll2007t iMuchovej et al.ll2010h . How- 
ever, there are significantly fewer sources in these regions 
compared to the cluster centers (a factor of ~ 10), and 
the contamination from these sources will be much less 
degenerate with the cluster SZ signal due to the relatively 
large spatial separation. Therefore, we expect the con- 
tamination from these sources to be significantly smaller 
than the contamination from the central radio galaxies. 
Consequently, we examine the effects of these sources in 
a statistical sense. 

First, the factor of 3 increase in the number of ra- 
dio sources will increase the corresponding noise fluc- 
tuations from undetected sources by a factor of ~ \/3. 
However, even with this enhancement the noise fluctu- 
ations from radio sources will still be sub-dominant to 
the noise fluctuations from dusty star-forming galaxies 
and negligible compared to the noise RMS in our images 
(< 0.1 mJy compared to ~ 1 mJy). Furthermore, we 
note that the increased number of radio sources in the 



direction of the cluster will systematically fill in the SZ 
decrement, which is referenced to the average signal level 
outside of the cluster, and bias our results. Based on a 
simple power-law extrapolation of the measured num- 
ber counts at ~ 140 GHz[3 which constrain the number 
counts above ~ 10 mJy, we estimate that the SZ decre- 
ment will be filled in by ~ 0.01 mJy/arcmin 2 due to 
the higher th an average number of radio sources t owards 
the cluster (IVieira et al.l 120101: iMarriage et~aTl l2011bl: 
\Planck early results XIH||2011D . This surface brightness 
is almost 3 orders of magnitude below the average SZ 
surface brightness towards the clusters in our sample and 
therefore should produce a ~ 0.1% bias in the average SZ 
signal recovered. Consequently, contamination from ra- 
dio galaxies outside of 1 arcmin in radius from the cluster 
center should generally be negligible. 

4. POSSIBLE CONTAMINATION FROM DUSTY GALAXIES 

In addition to radio emission, some of the cluster- 
member galaxies arc likely to have thermal dust emis- 
sion. Studies of clusters out to z — 0.8 in the mid- 
IR with Spitzer have shown that there can be signif- 
icant enhancements in the number of 24 /xm-detected 
sourc es within narrow redshift slices cente red on clus- 
ters (diEiFil] [2003; [Miicniiciril][2003). However, 
the total number of 24 yum objects detected towards 
clusters is < 10% above the number detected towards 
blank fields at similar sensitivities, due to the fact that 
the so urces are detected over a wide r ange of red- 
shifts ((GeacEetan [HI |FinEiEin[lQi3). Compared 
to 24 yitm observations, our data are significantly more 
sensitive to high-rcdshift galaxies and less sensitive to 
dust e mission from cluster-m ember galaxies (jBlain et al.l 
I2002t IBethermin et al.l 1201 ID , and so the effective en- 
hancement compared to a blank field will be even smaller 
in our Bolocam images. Furthermore, Herschel has re- 
cently published results characterizing the dust emission 
from 68 cluster BC Gs at frequencies as low as 600 GHz 
(|Rawle et al.H2012D . 15 of the BCGs arc detected in one 
or more Herschel bands, and grey-body fits of the Her- 
schel data result in extrapolated flux densities at our 
observing band of < 0.1 mJy for all but two of these 
sources. Based on our analysis of the contamination from 
central radio galaxies, even these brightest dusty sources 
would result in < 1% fractional changes in the SZ S/N. In 
summary, the dust emission sourced by cluster-member 
galaxies in our images is likely to be negligible, and we 
therefore do not account for it in this analysis. 

5. RADIO SOURCES DETECTED BY BOLOCAM OUTSIDE 

CLUSTER CENTERS 

The number of cluster member radio ga laxies drops 
quickly with radius from th e cluster center ([Coble et al.l 
120071 iMuchovei et al.ll2010D . and the blank sky descrip- 
tion of the source population becomes approximately cor- 
rect at these larger radii. Over the full set of 45 clusters 
we map a fairly significant area (~ 2.5 deg 2 ), which is 
enough to contain a handful of very brigh t radio sources 
(|Vieira et al.ll2010t IMarriage et al.ll2011bO . Fortunately, 
we are able to easily detect and subtract these sources 

16 Measurements at a wide range of frequencies indicate that 
a simple power law is a good approximation of the nu mber counts 
on th e low-flux side of the break in the counts (e.g., ITucci et al.l 

mm ). 
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TABLE 4 

Point Sources Detected by Bolocam 



cluster 


coordinates 
J2000 


NVSS 
1.4 GHz 


OVRO/BIMA 
28.5 GHz 


SZA 
30.9 GHz 


Bolocam 
140 GHz 


Abell 267 
Abell 383 
Abell 963 
Abell 2261 

MACS J2214.9 


01:52:54.6, +01:02:08.2 
02:48:22.1, -03:34:30.5 
10:17:14.2, +39:01:24.0 
17:22:16.9, +32:09:10.4 
17:22:23.8, +32:01:26.4 
22:14:39.3, -14:00:58.5 


4.6 ± 0.5 
54.9 ± 2.4 
1392.2 + 41.8 
24.3 ± 1.6 
126.7+ 4.4 
58.4+ 1.8 


7.55 + 0.24 

10.48 + 0.16 
107.10 + 0.58 


7.5 + 0.3 


9.8 + 2.3 

14.3 ± 1.7 
26.8 + 2.8 

8.7 + 1.1 
9.3 + 3.4 

51.4 + 5.2 



Note. — Point sources detected in Bolocam cluster observations with S/N > 4 in SZ-model-subtracted beam-smoothed images. The 
uncertainties given above fully account for degeneracies between the point source and the SZ model, and are therefore larger than the noise 
RMS in the beam-smoothed images. As a result, our constraint on the flux density of one source has a S/N < 4. Note that, in contrast to 
the central point sources selected from NVSS, 5/6 of these sources appear to be flat-spectrum sources. This is due to the fact that most 
stccp-spcctrum sources are too dim to be detected in our 140 GHz data. 



from our data because they have little or no degener- 
acy with our template of the SZ signal. We search for 
these sources by selecting pixel excursions with a S/N 
> 4 in maps that have had the best-fit SZ template 
subtracted and are beam-smoothed for improved point 
source extraction. We find 15 such excursions, 12 with 
positive flux density and three with negative flux den- 
sity. Six of these positive flux candidates have counter- 
parts in NVSS, and we estimate the 140 GHz flux density 
of each of these sources by fitting a point source tem- 
plate to our maps us ing the NVSS coordinates as a prior 
(|Condon et al.l [l998h . The dimmest of the sources with 
NVSS associations have flux densities of < 10 mJy in 
our images, providing a rough estimate of our detection 
threshold (see Table |4}. After removing the six sources 
with NVSS associations, we are left with nine unasso- 
ciated pixels (six positive and three negative). All of 
these unassociatcd pixels have S/N just above four. Sta- 
tistically, based on Gaussian noise, we only expect ~ 1 
noise excursion with S/N larger than four, and it is un- 
clear why we have nine such excursions. However, we do 
note that most of our images have significant SZ signal 
in them (S/N > 10), and relatively modest differences 
between our SZ model and the underlying cluster profile 
could therefore combine with noise excursions in some 
clusters and cause such high S/N pixels in our best-fit- 
SZ-model-subtracted images. In addition, the noise in 
our images is not perfectly white, and large-angular-scale 
noise could also be responsible for some of these S/N > 4 
pixels. We do not attempt to further account for these 
nine unassociated pixels with S/N > 4. 

6. CONCLUSIONS 

Using a combination of 140 GHz Bolocam data, 
1.4 GHz NVSS data, and ~ 30 GHz data from 
OVRO/BIMA and SZA, we have studied the 140 GHz 
emission from cluster-member galaxies in a sample of 45 
massive clusters. In agreement with previous results, we 
find that cool-core clusters preferentially host more radio 
galaxies and that all of the cluster-member galaxies have 
steep spectra. On average, we find that the spectral in- 
dex between 30—140 GHz is consistent with, but slightly 
steeper than, the spectral index between 1.4—30 GHz. As 
a result, we find a multiplicative factor of 0.85 ± 0.16 be- 
tween the sample-averaged 140 GHz measured flux den- 
sities and the predicted flux densities based on an ex- 
trapolation from the lower frequency data. We further 
find that galaxy to galaxy variations lead to a ~ 30% 



scatter in 140 GHz flux densities extrapolated from data 
below 30 GHz. Based on an extrapolation using our best- 
fit steepening of the spectral index, which results in a 
15% correction, we find that only ~ 1/4 of the clusters 
contain enough radio emission to produce more than a 
1% bias in the optimally filtered peak SZ signal from 
the cluster, and the maximum contamination in a single 
cluster is ~ 20%. However, we do note that the con- 
tamination in cool-core systems is significantly enhanced 
compared to the sample average. These results roughly 
match expectation s from lower-frequency measurements 
and simulations dLin et al.l 120091: iSehgal et all 120101: 
lAndersson et al.l 120 111) and indicate that the amount of 
radio contamination in high mass clusters is small com- 
pared to achieved noise levels in SZ measurements ob- 
tained at 140 GHz. 
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